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ABSTRACT 

We propose the existence of ultracompact minihalos as a new type of massive compact halo object 
(MACHO) and suggest an observational test to discover them. These new MACHOs are a powerful 
probe into the nature of dark matter and physics in the high energy Universe. 

Non-Gaussian energy-density fluctuations produced at phase transitions (e.g., QCD) or by features 
in the inflation potential can trigger primordial black hole (PBH) formation if their amplitudes are 
S £ 30%. We show that a PBH accumulates over time a sufficiently massive and compact minihalo to 
be able to modify or dominate its microlensing magnification light curve. Perturbations of amplitude 
0.03% ;$ S ,$ 30% are too small to form PBHs, but can nonetheless seed the growth of ultracompact 
minihalos. Thus, the likelihood of ultracompact minihalos as MACHOs is greater than that of PBHs. 
In addition, depending on their mass, they may be sites of formation of the first Population III stars. 

Ultracompact minihalos and PBHs produce a microlensing light curve that can be distinguished from 
that of a "point-like" object if high-quality photometric data are taken for a sufficiently long time 
after the peak of the magnification event. This enables them to be detected below the stellar-lensing 
"background" toward both the Magellanic Clouds and the Galactic bulge. 
Subject headings: early universe — dark matter — gravitational lensing — Galaxy:halo 



1. INTRODUCTION 

The nature of dark matter is one of the fundamen- 
tal unsolved questions in Cosmology. CMB anisotropy 
data indicate that about 20% of the Universe is com- 
posed of non-baryonic dark matter and 4% is in baryons. 
The most popular dark matter candidates are weakly in- 
teractive massive particles (WIMPs) and a great effort 
is under way to detect them. Ongoing experiments us- 
ing ground based dark matter detectors are significantly 
restricting the paramete r space for some f avored super- 
symm etric models (e.g.. lAngle et al.l 2008; A hmed et al.l 
120091) . 

A less popular - but physically well-motivated - dark 
matter candidate is primordial black holes (PBHs). 
These are black holes predicted to form during the ra- 
diation dominated era if the Universe had regions with 
energy-density fluctuations of amplitude £ 30%. PBHs 
can form with a range of masses that span many decades, 
from the Planck mass (10~ 5 g) to thousands of solar 
masses. 

A fully relativistic approach is necessary to model the 
formation of PBHs. However, the Newtonian approxi- 
mation elucidates the basic concept of why PBHs can 
form before large scale structures appeared in the Uni- 
verse. During the radiation era the cosmic Jeans mass 
approaches the Horizon mass, which is also of the same 
order of the mass of a black hole with density equal to 
the mean cosmic value. This means that before matter- 
radiation equality relatively small perturbations on Hori- 
zon scales may become self gravitating and they begin the 
collapse with an initial density that is already very close 
to the black hole regime. Thus, the mass of PBHs de- 
pends on the time of their formation because it is of the 
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order of the mass of the particle Horizon at that time. 

However, fluctuations of amplitude of about 30% - re- 
quired for PBH formation - are very large when com- 
pared to the r.m.s amplitude of perturbations from in- 
flation (about 0.001%) that lead to the formation of nor- 
mal galaxies. Thus, such large perturbations would be 
exceedingly rare, unless they can be produced by physi- 
cal processes taking place during phase transitions (e.g., 
topological defects, bubble nucleation, softening of the 
cosmic equation of state, etc) or by ad hoc features in 
the inflation potential. Due to theoretical uncertainties 
in the physics of the high-energy Universe we do not know 
the typical amplitude of perturbations created at phase 
transitions. Hence, theoretical modeling cannot tell us 
whether PBHs exist, whether they are the bulk of dark 
matter or whether they are only an academic curiosity. 
We must rely on observations to probe them as dark 
matter candidates. 

Even if PBHs are only a fraction (e.g., < 10%) of 
the dark matter, during matter domination, they can 
grow by up to two orders of magnitude in m ass through 
the a cquisition of large dark matter halos (jMack et al.l 
120071 hereafter MOR07). The dark halo is instrumen- 
tal in increasing the ability of the PBHs to accrete gas, 
in boosting their X-ray emission, a nd thus in modify- 
ing the cosmic ionization history (Ric otti et al.l 120081 . 
hereafter ROM08). In this paper we consider whether 
the dark halo is sufficiently massive and compact to 
produce observable modifications of microlensing light 
curves of PBHs or even dominate the inferred mass of 
the lens. Most importantly, we propose a novel type of 
non-baryonic MACHO. The idea is a corollary of models 
of PBH formation: if the amplitude of the fluctuations 
that may induce collapse of PBHs is smaller than the 
required threshold of ~ 30%, then even though the PBH 
is not formed, the perturbation will nonetheless seed the 
growth of ultracompact minihalos observable using mi- 
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crolensing experiments. Due to the lower overdensity 
threshold required to form ultracompact minihalos, their 
existence is several orders of magnitude more likely than 
that of PBHs. 

Observational evidence for (or constraints against) 
PBHs is often soug ht in conflicting clai ms for (MACHO: 
Alcock et ail |2000fi or against (EROS : iTisserand et all 
2007t OGLE: IWvrzvkowski et al.ll2009ft detection of mas- 
sive compact halo objects (MACHOs) in microlensing ex- 
perime nts toward the Ma gellanic Clouds (MCs) . For ex- 
ample, lAlcock et al.ll2000l originally claimed that roughly 
20% of Galactic dark matter was in the form of MA- 
CHOs of characteristic mass 0.4 M@, which made PBHs 
a good candidate (since stars of this mass would easily be 
seen). However, ROM08 argued that such a population 
of PBHs would emit X-rays and produce distortion of 
the CMB spectrum that are incompatible with COBE- 
FIRAS data. 

Here, however, we adopt a very different orientation. 
First, we reopen the possibility that PBHs can be iden- 
tified with microlenses detected toward the MCs by not- 
ing that the PBH contribution to the microlens mass is 
small compared to that of accreted dark matter. Hence 
the X-ray signature is likewise small. Second, we argue 
that these PBH-seeded compact objects can have most of 
their mass outside the lensing Einstein radius, meaning 
that the objects might be cosmologically very important 
even if naive interpretations of the lensing results con- 
strain their density to be < 10% of Galactic dark matter 
(which would be consistent with both experiments claim- 
ing upper limits). Third, and most important, we argue 
that because they are embedded in accreted dark mat- 
ter, PBHs generate a lensing signal that is qualitatively 
different from that due to stars and other "point-like" 
objects. Hence, PBHs would be detectable (and distin- 
guishable from stars) in microlensing experiments even 
if their lensing rate was equal to or below that of stars. 
Moreover, even if the amplitude of early-universe per- 
turbations is too small to create PBHs, it may still be 
large enough to generate compact minihalos that would 
give rise to lensing signatures. Like the PBH-minihalo 
signatures, these could also be robustly distinguished 
from garden-variety microlensing due to stars. Thus, we 
propose an observational test for PBHs and other early- 
universe perturbations that is far more sensitive and ro- 
bust than any previously contemplated. 

The plan of the paper is as follows. In § [2] we dis- 
cuss the formation and derive the density profile of ul- 
tracompact minihalos. In § [3] we calculate the magni- 
fication light curve of microlensing events due to PBHs 
embedded in their minihalos and ultracompact miniha- 
los without PBHs. In § 0] we discuss the prospects of 
detecting PBHs and ultracompact minihalos below the 
stellar-lensing "background" toward both the Magellanic 
Clouds and the Galactic bulge. Finally, in §[5] we present 
the discussion and conclusions. 

2. DENSITY PROFILE OF ULTRACOMPACT MINIHALOS 

It is generally accepted that the mass of PBHs 
does not increase sig n ifican t l y after their formatio n 
dZerDovich fc Novikovi fl967t ICarr fc Hawking) Il97l . 
However, any locally overdense region in an expanding 
universe seeds the formation of dark matter structures. 
PBHs are local overdensities in the dark matter distri- 



bution (made either of WIMPs or smaller mass PBHs), 
hence they seed the growth of spherical halos. 

The theory of spherical gravitational collapse in an 
expanding universe (i.e., assuming radial infall), als o 
known as secondary infall theory ( Bcrtschi ngerl [1985T ) . 
predicts that during the matter-dominated era, the dark 
halo grows as i 2 / 3 oc (1 + z)~ l . During the radiation- 
dominated era the halo growth is of order unity, thus the 
halo mass grows with redshift as 

M h {z) = M pbh (^^-) , (1) 

where z eq 3500 is the redshift of matter-radiation 
equality (MOR07). After z ~ 30, the growth of the 
dark minihalo depends on the environment. If the PBH 
evolves in isolation it can continue to grow, otherwise it 
will either stop growing or it will lose mass due to tidal 
interactions as it is incorporated into a larger galactic 
halo. 

2.1. Revisiting Secondary Infall Seeded by PBHs 

In previous work on halo growth seeded by PBHs, 
MOR07 have assumed a point-like overdensity (i.e., the 
PBH) accreting gas and dark matter from a uniform den- 
sity expanding universe with p = Q m p C rit (the growth 
of the halo mass during radiation epoch is negligible). 
However, this assumption is correct only if the mass of 
the perturbation that originates the PBH is equal to 
the PBH mass. This is typically incorrect. We expect 
Mpbh ~ Mhov, only if the perturbation that creates the 
PBH has amplitude much larger than the critical ampli- 
tude for collapse. 

One dimensional fully-relativistic simulations of grav- 
itation collapse of perturbations on Horizon scales show 
that the critical overdensity needed to trigger PBH col- 
lapse is » w , wher e P = wpc 2 is the cosmic equation 
of state (|Carrlll975t iGreen et al1l2004| ). During the ra- 
diation era w — 1/3 ~ 30%. The mass of PBHs is typi- 
cally smaller than the mass within the particle Horizon at 
the redshift of its formation, zf. M p },h = fHorMjior(zf), 
where fuor < 1 depends on the amplitude of the per- 
turbation relative to the critical value w. Perturba- 
tions with amplitude significantly larger than w produce 
ftior ~ 1) while perturbation that are near the criti- 
cal value produce PBHs with masses much smaller than 
Mhot- In addition, the total mass of the perturbation 
can be 5m > Mn or (zf) if it has wings of subcritical am- 
plitude extending outside of the Horizon at Zf. Partially 
for t his reason it has been suggested by several authors 
(e.g.. lDokuchaev et aLlfeOO^ ICarrl[2005l IChisholmll2006| ) 
that PBHs have high probability of forming clusters or 
binaries. 

According to the theory of secondary infall any mass 
excess 5m > within a sphere of radius R seeds the 
growth of a minihalo. Here 5m = M(R) — 4ir/3£l m p cr R 3 
includes the mass of the PBH and any overdensity 
surrounding it. We ran several simulations (see Sec- 
tion [2TTT]) in which the PBH is located in an overdense 
region confirming that the mass of the minihalo is pro- 
portional to the mass of the overdense region 5m: 

M h {z) = 5m(±±f) . (2) 
V 1 + z eq J 
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If the accretion has spherical symmetry (i.e., radial in- 
fall), the dark halo develops a self-si milar power-law den- 
sity profile p oc R~ v with v ~ 2.25 (|Bertschingerlll985ft 
truncated at a halo radius 



R h = 0.019 pc 



M h (z) 



1 Mi 



^0 



1/3 



l + Z 

1000 



(3) 



where Rh is a bout one third of the turn-around radius 
(|Ricottil 120071 ROM08). Hence, the cumulative mass of 
the halo is M h {R) = M h (z)(R/R h (z)) 3/i . 

2.1.1. Simulations 

We have run ID simulations of secondary accretion 
that confirm the results expected from the theory of sec- 
ondary infall. The code we used is a modified version of 
the code described in MOR07. The modification allows 
for shell crossing and crossing of the singularity at R = 0. 

We ran simulations with different shape of the overden- 
sity (top-hat, triangular, etc), different overdensity, and 
different size. After a transient phase the density profile 
always settles to a power law with slope 2.25. The mass 
profile of the halo depends on 5m as in equation ([2|). 
However, if the initial amplitude of the perturbation at 
equality is too small for the perturbation to grow to 5 ~ 1 
by a given redshift, the collapse has 2 phases: during the 
first phase the halo is less dense and in the second phase, 
after i5 ~ 1 the density increases to its final value. We es- 
timate the critical amplitude required for the formation 
of an ultracompact minihalo as the one that leads to the 
second and final phase of the collapse at z ~ 1000. This 
is because halos that reach the final density profile at 
lower redshifts are likely to have a shallower inner core 
due to the larger angular momentum of accreted dark 
matter leading to some degree of non-radial infall. 

2.1.2. Lensing 

The microlensing magnification by the minihalo and 
PBH depends on the mass profile projected on a plane. 
In Section [3] we find that Mj_(R) = l.350M h (R) for a 
density profile with v = 2.25. Hence, from equations |1|- 
©: 



M X (R) = 0.058 M Q 



5m 

1 Mr. 



3/4 



R 



AU 



3/4 



(4) 



The reason that we parameterize the radius in equa- 
tion Q in units of 8 AU is that the Einstein ra- 
dius of a lens with mass Mi ens is about i?E ~ 
8 AU(M Ie „ s /l M e ) 1/2 . The mass of the lens is M iens « 
M ± (R E ) + M pbh . If M ± {R E ) < M pbh , then M lens « 
M pbh and 



M±(R E ) 



M, 



pbh 



5m 
1 M G 



3/4 



1 Mr, 



-5/8 



(5) 



In Section [3] we will show that if M±(R E ) Z 30%M pbh 
the difference of the light curve of the lensing event from 
that of a point mass with M; eras = M pb h is measurable. 
If 5m — M pb h the effect of the minihalo in modifying 
the magnification light curve of the PBH is probably 
not observable assuming realistic measurement errors. 
However, we have already pointed out that typically 
5m ~ M Hor > Mpbh- If 6m /Mpbh > 9 (or f Hor < 0.11 



assuming 5m = MHor) the signature of the minihalo is 
observable in the lensing light curve. 

In the limit of Mj_(R E ) > M pbh , then M lens m 
M±(R E ) and we find 



Mi 



lens 



1 Mr 



5m 



44.5 Mr 



6/5 



(6) 



In this case the mass of the lens is independent of the 
PBH mass. In the next section we consider the case in 
which the PBH does not form: M pb h = 0. This case is 
particularly interesting because the formation of ultra- 
compact minihalos is statistically more likely than the 
formation of PBHs. 

2.2. Ultracompact Minihalos without PBHs 

Dark matter density fluctuations that seed the forma- 
tion of large scale structure and galaxies, have an am- 
plitude 5 ~ 10~ 5 when they enter the Horizon, roughly 
independently of their mass (for scale invariant pertur- 
bations). Here we focus on small mass fluctuations that 
enter the Horizon during the radiation era. The rela- 
tivistic component of the (adiabatic) perturbations stops 
growing rather quickly after it enters the Horizon be- 
cause its mass becomes smaller than the Jeans mass 
(Mj(z) ~ Mhot{z) for z > z eq ). The dark matter 
component grows only logarithmically because of the 
Meszaros effect. Overall, small mass linear perturbations 
grow by 2-3 orders of magnitude from the time they en- 
ter the Horizon to matter-radiation equality, depending 
on their mass. A relatively rapid initial growth when the 
perturbation enters the Horizon is followed by a slow log- 
arithmic growth. For the mass range of interest to us, the 
Horizon mass increases by 10 or more orders of magni- 
tude from the time the perturbation enters the Horizon 
to to the redshift of equality (MH or (z eq ) ~ 10 15 M ). 
For example, scale invariant density fluctuations with 
mass typical of dwarf galaxies (10 8 — 10 9 M Q ) grow to 
an amplitude 5 ~ 10 -2 by redshift of matter-radiation 
equality. They collapse when 5 ~ 1.68 at z V i r ~ 15 fol- 
lowing further growth by a factor z eq /z V i r ~ 200 during 
matter-domination era. 

Dark matter perturbations that enter the Horizon with 
amplitude £ 0.1% - i.e., ~ 100 times larger than the ap- 
proximately scale invariant perturbations from inflation 
- collapse before z — 1000. These perturbations are too 
small to form PBHs but can seed the growth of ultracom- 
pact minihalos that accrete from a nearly uniform Uni- 
verse, before the epoch of formation of the first galaxies 
at z ~ 30. Thus, the angular momentum of the accreted 
dark matter is small and the infall quasi-radial. 

The aforementioned arguments lead us to conclude 
that the formation of ultracompact minihalos is more 
likely than the formation of PBHs. The threshold over- 
density required for their formation at z ^ 1000 is ~ 
0.1%, compared to « 30% required for PBH formation. 

Similarly to the case for PBHs, observational con- 
straints or the discovery of ultracompact minihalos is a 
powerful probe of high-energy physics in the early Uni- 
verse. Microlensing experiments can help constrain the 
amplitude of perturbations produced during two recent 
phase transitions: the QCD (quark-hadron) phase tran- 
sition with Mhot ~ 1 M Q and the e + — e~ annihila- 
tion epoch with Mjjor ~ 10 5 M Q . Perturbations with 
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Fig. 1. — (Left) Time sequence of the growth of the mass profile of a spherical halo due to secondary infall. The halo growth is seeded by 
a density perturbation with overdensity S(z eq ) = 10 and excess mass 5m = 100 Mq. The solid lines show the time evolution of the mass 
profile from z = 1000 to z = 5. The long dashed line delineate the lens equation. The two dotted lines show power law mass profiles with 
log slope 0.75 producing gravitational lenses with Af; ens = 0.01, 1 Mq. (Right) Same as the left figure for a perturbation S(z eq ) = 0.1 and 
5m = 100 M . 

o'dm.Rh < (GMiensRE) 1 / 2 ■ This inequality is satisfied 
for 



masses in this range can be constrained by microlensing 
experiments because the mass of the lens, given by equa- 
tion ©, is in the range probed by current microlensing 
experiments. 

2.3. Angular Momentum of Accreted Dark Matter 

The angular momentum of accreted dark matter de- 
termines the density profile of the dark halo enveloping 
a PBH and whether the mass of a PBH can grow sub- 
stantially by accreting a fraction of its enveloping dark 
halo. ROM09 found that angular momentum of dark 
matter accumulated by PBHs with masses M p bh > 1000 
is so small that a fraction of the dark matter is directly 
accreted by the PBH. But for smaller mass PBHs and 
ultracompact minihalos the angular momentum is larger 
and the accretion can become non radial in the inner 
parts of the halo. The v = 2.25 power law profile of the 
density of minihalos results from the assumption of radial 
infall. If the accreted material has some angular momen- 
tum the density profile will be shallower near the center 
and will contain less mass. Here we check whether the 
assumption of radial infall is justified in the inner parts 
of the halo profile enclosing a mass comparable to the 
lensing mass (i.e., within the Einstein radius). 

The angular momentum of the dark matter accret- 
ing onto PBHs can be estimated from equation (19) in 
ROM09, which gives the mean values of the velocity dis- 
persion within a comoving volume of radius equal to the 
halo turnaround radius. 



5m < 1350Mp 



© 



-i . \ 3.46 / , , r 

1 + z \ l Mi, 



1000 



1 Mr. 



(8) 



In equation ((8|) there is a steep dependence on redshift. 
A value of redshift z ~ 1000 is most appropriate for 
probing the inner parts of the density profile because the 
inner m ass is domina t ed by material accreted at high 
redshift. iBertsching er (1985) has shown that the growth 
of the halo is self similar because the material accreted at 
late times from outer shells does not change the density 
profile in the inner parts (despite the fact that these carry 
most of the halo mass). This is because the outer shells 
have high speed when they reach the halo center and 
spend little time there. Thus, the mass profile near the 
center of the minihalo - within the Einstein radius - is 
built by material accreted at redshifts z ^ 1000. Let us 
now rewrite equation ([8]) for two special cases: 

1) If Mi ens w M pbh , assuming 5m = M pbh / f Hor we 
find that the infall is quasi-radial at Re if 



7.4 x 10" 



, . N -3.46 / u \ -0.22 
1 + Z \ I M p bh 



1000 



1 M m 



< fHor < 1. 



(9) 



l + z 
1000 



Mh 

1 Mr. 



0.28 



(7) 



2) If Mi ens ~ M± from equation © we find that the 
infall is quasi-radial at Re if 

-7.4 

5m > 2.9 x 10~ 2 M ' ' 



l + z 



1000 



(10) 



with (Jdmfi = 1.4 x 10~ 4 km s~ . 

Applying conservation of angular momentum we find 
that the rotational (i.e., tangential) velocity of the gas 
at a distance r from the black hole is v(r)r — o~dmRh, 
where Rh is the halo radius. If the velocity is smaller 
than the Keplerian velocity in the proximity of the Ein- 
stein radius, then the accretion is quasi-spherical; if vice 
versa, a shallower core can form. The Keplerian ve- 
locity at radius R is {GM{R)/ 'R) 1/2 . Thus, the ac- 
cretion is quasi-spherical within the Einstein radius if 



2.4. 



Profile Steepening due to Adiabatic Compression 
and Population III Formation 

In this section we discuss the role of baryons on shap- 
ing the density profile of ultracompact minihalos and 
whether the first Population III stars can form in their 
centers. The discussion here will be only qualitative as 
the simulations required to obtain quantitative results 
are beyond the scope of the present paper. 

Due to gas cooling and the concentrated dark mat- 
ter profile of ultracompact minihalos, baryons may sink 
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dissipatively by several orders of magnitude in radius to- 
ward the halo center. The infall perturbs the underlying 
dark matter distribution pulling it inward and creating 
an even smaller and denser core than would have evolved 
without baryon condensation. For the case of ultracom- 
pact minihalos, it is most appropriate to assume perfectly 
radial orbits (but the same equations are valid assuming 
only circular orbits). Since Mh(R) varies in a self-similar 
fashion, the quantity Rm ax M h {R max ) , where R max is the 
maximum radius of the radial orbit, is conserved during 
the "slow" contraction of the gas component. Defining 
the fraction of baryonic mass F = M^/Mh, the invariant 
of the dark matter particle orbits implies 

R[M b (R) + M h (R)} = R^'h{n) = ^iMhM (ii) 

1 — r 

which can be solved iteratively for the final dark mat- 
ter mass distribution M^R), given the initial total mass 
distribution Mh i(Rj) and the final bar yon mass distri- 
bution M b (R) (jBlumenthal et alJ Il986t ). The equality 
Mh.i = (1 — F)Mi(Ri), further assumes that initially 
F(Ri) = const. In Figure [2] we show the compres- 
sion of the dark matter profile due to baryonic con- 
traction obtained by solving iteratively equation (jTTJ) . 
The initial minihalo mass profile is a power law Mh,i = 
M h (z)(R/R h ) - 75 with M h (z) from equation (0) assum- 
ing 5m = 1 M Q and z = 10 (dotted line). The final 
dark matter profiles are shown by the solid and dashed 
lines for different configurations of the final baryon dis- 
tribution. We assume that a fraction of the baryons 
fbar(core) — 0.01, 0.001 settles to a constant density core 
with radii R CO re = 0.5 AU and 8 AU. The value of R cor e 
is not important as long as R core < Re- The mass of the 
lens is given by the intersection of the final dark mat- 
ter mass profile with the long dashed line, which traces 
the lens equation M lens /(1M Q ) « [R E / (8AU)} 2 . The 
mass of lens increases only if the gas mass within Re for 
the initial mass profile is comparable or dominates initial 
dark matter mass. If later, the baryons are removed from 
the center of the halo on a short time scale when com- 
pared to the orbital period of the dark matter particles 
(for instance due to reionization, SN explosions, etc) the 
cuspy dark matter profile may become less steep, but it 
will not go b ack to the original p r ofile before the baryon 
compression (IGnedin et al.ll200l ISellwood fc McGaughl 
120051) . 

The next step consists in determining whether a frac- 
tion of the baryons in the minihalo can become suffi- 
ciently concentrated to dominate the mass within Re at 
some point during the minihalo evolution. We will con- 
sider separately the case of PBHs embedded in minihalos 
and ultracompact minihalos without PBH. 
1) For PBHs accreting at high-z from the IGM, the gas 
mass within the Einstein radius of the lens is proportional 
to the dimensionless gas accretion rate fa = M g /MEd, 
where Me<i is the Eddington rate. For PBHs with mass 
;$ 500 M the gas accretion is quasi-radial. Thus, as- 
suming Mi ena <~ Mpbh, the gas approaches free fall ve- 
locity v ~ (G Mi ens / R) 1 / 2 - From mass conservation 
ATrp(R)R 2 v = mM E d, we find that 

MJR E ) « /Aliens \ 1/4 , s 

gK ' ~ 5.4 x 10 _8 m . 12 




It)- 8 10- 6 10"* 10" 3 0.01 0.1 1 
R (pc) 



Fig. 2. — Dark matter mass profile of minihalos before (dotted 
line) and after (solid and dashed lines) baryonic compression. The 
initial mass profile is for a minihalo at z = 10 with 5m = 1 Mq. We 
show the mass profile for a set of final baryon configurations. For 
simplicity here we assume that the baryons have a constant density 
core with radius R CO re = 0.5 AU (solid lines) and 8 AU (dashed 
lines) and with gas mass that is a small fraction ff, ar (core) = 
0.01, 0.001 of M h (see labels). 

Thus, in this regime, the gas mass never dominates the 
potential within Re- However, equation (| 1 2[) does not 
apply in the following cases: i) PBHs more massive than 
Mybh > 10 5 Mp), for which the gas is self-gravitating 
(|Ricotti|[2007l ); ii) PBHs with mass > 500 - 1000 M 
accreting gas at z < 100, for which the gas has sufficient 
angular momentum to form an accretion disk. Thermal 
feedback effects, which are important at z < 100, will 
tend to increase the value of the PBH critical mass for 
disk formation (ROM08). An accretion disk may also 
form around PBHs that are able to accrete gas efficiently 
from the ISM in galaxies. However, to sustain efficient 
gas accretion the relative velocity of PBHs with respect 
to the ISM must be subsonic. Only in this latter case, if 
the accretion rate is near Eddington, the accretion disk 
may dominate the mass profile within Re- We conclude 
that compression of dark minihalos around PBHs due to 
baryon dissipative collapse is in most cases negligible. 
2) The first stars are thought to form i n rare 3 a per- 
turba tions, of mass ~ 10 5 M a atz~ 30 (|Tegmark et alJ 
Il997| ). Ultracompact minihalos, because they are already 
in place at z ~ 1000, may form the first stars in smaller 
mass halos of 10 4 M Q at redshifts 200 < z < 1000. The 
simplest and naive condition for star formation is given 
by the requirement that the minihalos must have T V i r > 
T cm b- the gas would not be able to condense if it were 
initially colder than the tem perature of the CMB (see 
Fig 6 in Teg mark et al.llT997t ) . Thus, ultracompact mini- 
halos with 5m > 100 M Q are able to host the formation 
of Population III stars, and the baryon dissipation can 
increase the concentration of the dark matter halo. 

However, the naive estimate is likely too conservative 
because ultracompact minihalos have a power law density 
profile with v — 2.25, thus are much more dense than nor- 
mal galaxies (e.g., assuming NFW profile). Their steep 
mass profile leads to significantly larger gas densities in 
their cores. Assuming an isothermal equation of state 
with T g = Tcmb (a good approximation at z > 100 due 
to Compton cooling/heating), the gas density profile in 
hydrostatic equilibrium in the potential of an ultracom- 
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pact minihalo is 

Pb(r) =p h exp^ a 



R_ 

Rh 



-1/4 



1 



(13) 



where a = Ar b /r h (|Ricottill2009h . For T g = T C mb we get 
a = (M h /240 M Q ) 2 / 3 « 0.054<5m 2 / 3 (l + z/1000)" 2 / 3 . If 



we estimate pb{R) at R = Re, where 



^ = 1.6x10-% lf 
Rh \IM Q 



13/15 



1 



1000 



4/3 



we find that the core gas density is 



Pb(Ri 



Pb ex P S a 



R_ 

Rh 



-1/4* 



p b io°' 37( " 



(14) 



> \ 1000 > 



(15) 



The gas density in the core becomes increasingly large 
with decreasing redshift. Equation (fT5|) is a good esti- 
mate of the gas density in the core as long as the gas 
temperature is close to or lower than the CMB temper- 
ature (i.e., z > 10 — 30) and pb(Rs) < Pdm(RE)- 

If the gas in the core is able to form H2 and cools in 
less than a Hubble time, it may collapse and form the 
first Population III stars. The large initial density in 
the core increases the ra te of formation of H2 with re- 
spect to the estimates by iTegmark et al.l (|1997ft for nor- 
mal galaxies (although the reduced fractional ionization 
and H~ abundance in the high density core have the 
opposite effect, these do not dominate). Thus, ultracom- 
pact minihalos of mass significantly smaller than 10 4 M Q 
at z ~ 30 may be able to form dense gas cores or Popula- 
tion III stars. It is therefore possible that minihalos with 
5m 1 Mq, as the one shown in Figure^ may be com- 
pressed by baryon collapse and produce lenses with much 
larger masses than given by equation ([6]). Quantitative 
calculations of the gas collapse require the integration of 
a set of time-dependent equations of the chemical reac- 
tion network for H2, in the evolving gravitational poten- 
tial of the minihalo. These calculations are beyond the 
scope of the present paper. 

3. MICROLENSING LIGHT CURVE 

To calculate the magnification of minihalo MACHOs 
(initially without PBHs at their center), we consider den- 
sity profiles of the form p = Cr~ v , where C is set so that 
the projected mass M± within a radius Re is just equal 
to mass required to form an Einstein ring at R E . That 
is, 



M ± {R)= / 2ndbb 
Jo 

AGM ± {R E ) 



dzC(b 2 + z 2 Y 



(16) 



R p — 



D 



^-k-k) <17) 

where is the distance from Earth to the minihalo lens 
and Ds is the distance to the microlensed source. Note 
that M±(R) is related to the mass inside a sphere of 
radius R by 



M ± (R) (l/2)![(^-3)/2]! 
M(R) \{v-2)/2}\ 

which is 1.350 for v = 2.25. Hence, M ± (R) oc R s ~ u 



(18) 



We then use Newton's method to solve the lens equa- 
tion, which gives the image position 9i as an implicit 
function of the source position 9s, 



?S = 



4GM ± (D L 9i 



1 



1 



7^-7H- (19) 



Newton's method automatically returns the radial 
derivative of M± , and so allows calculation of the mag- 
nification, 



A = A A 



A+ = 



89i di 
09 s 9 S 



(20) 



For the case that the minihalo is seeded by a PBH, we 
1 simply add a point mass to the center of the halo profile 
'and repeat the above procedure. 

We construct simulated light curves (with errors pro- 
portional to the square root of the magnification) and 
proceed to fit the simulated flux F to a standard 5- 
parameter point-lens microlensing model of the form 



A(u) 



F = f s A(t;t ,u ,t E ) + f b 
u 2 + 2 



u\lu 2 +4' 



«(*) = 




where t E is the Einstein radius crossing time, uo is the 
impact parameter in units of the Einstein radius, to is 
the time of maximum (closest approach) , f s is the model 
source flux, and /& is background light that lies in the 
photometric aperture but does not participate in the 
event. We adopt the conservative approach of allow- 
ing negative blending (which would be unphysical) for 
two reasons. First, the great majority of microlensing 
events are blended because they are in crowded fields. 
This blended light typically absorbs any negative blend- 
ing from a spurious fit, rendering the latter undetectable. 
Second, even truly isolated sources can suffer negative 
blending because of poor determination of the sky. 

The two examples in Figure 3 show deviations from 
standard microlensing of somewhat different degree. 
Whether either of these would be detectable in practice 
would depend on the quality of the experiment. However, 
the crucial point is that they are detectable just based on 
their structure. When microlensing dark matter searches 
were proposed by Paczynksi (1986), their major selling 
point was that if the dark matter were composed of MA- 
CHOs, then microlensing events would be 25 times more 
common than if it were not. Hence the event rate itself 
would unambiguously settle the question. Then, when 
MACHO reported an event rate that was the root-mean- 
square of the prediction for MACHOs and the prediction 
for stars, the result became completely ambiguous. Re- 
cent work by EROS and OGLE has seemed to show that 
the event rate is compatible with that expected from 
stars. If confirmed, this would have seemed to be the 
end-game for microlensing dark-matter searches: even if 
there were dark-matter objects, they would be too few to 
detect as an excess event rate (without enormously larger 
experiments) and, in any event, would be cosmologically 
uninteresting because their density was so far below the 
dark-matter density. 

But Figure 3 shows that minihalo MACHOs have light 
curve signatures that are different from those of stars and 
therefore can be detected even if they are only as common 
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Fig. 3. — (Left). Magnification light curve (top panel) of a PBH of 0.5 Mq including its minihalo with mass within = 8AU of 1 Mq. 
Here the source impact parameter and crossing time in units of the Einstein radius are 0.05 and 30 days, respectively. The bottom panel 
shows the residual with respect to the magnification curve of a point mass. The synthetic data were taken for 10 Einstein ring crossing 
times (about 10 months in this example). (Right ). Same as the left figure for an ultracompact minihalo (without PBH) with M; ens = 1 Mq. 



as stars, or indeed, even if they are far less common. 
Moreover, as we have discussed in Section [U minihalo 
MACHOs could make up all the dark matter and yet 
have a microlensing optical depth that is one or several 
orders of magnitude below that expected for point-lens 
dark matter. 

We note that there are several observed microlensing 
effects that produce deviations from the standard mi- 
crolensing form given by equation (|22[) . For example, 
parallax effects due to the Earth's acceleration, xallarap 
effects due to source acceleration from an orbiting com- 
panion, light curve caustics due to binary or planetary 
companions, and finally finite-source effects that occur 
if the lens transits the source. However, except for the 
last, all of these are generically asymmetric, whereas the 
light curves shown in Figure 3 are symmetric. Moreover, 
although finite-source effects are symmetric, they have a 
very specific form which is unlike the form due to mini- 
halo MACHOs. Furthermore, they are almost completely 
restricted to extremely high-magnification events, which 
are rare. 

Now it is true that the other effects could, in particular 
cases, be very symmetric but this is highly improbable, 
and furthermore these effects, even in their more typi- 
cal asymmetric form, are not all that common. Hence, 
the level of contamination in minihalo-MACHO searches 
would be extremely low. 

4. MACHO SEARCHES TOWARD THE BULGE 



Originally, iPaczvnskil ([1986) proposed MACHO 
searches toward the LMC not because this line of sight 
had the highest expected number of MACHO microlens- 
ing events, but because it had very little contamination 
from non-MACHO events. However, if the MACHO and 
non-MACHO events can be distinguished based on their 
light curves, one should really focus MACHO searches on 
the line with the highest rate of dark-halo events and the 
highest-quality light curves (needed to decisively charac- 
terize the deviations from point-lens microlensing). This 
is the Galactic bulge, where currently about 800 events 



are discovered per year and where the sources are 6 times 
closer, so 36 times brighter (for similar luminosities). It 
is true that the fraction of these due to MACHOs is 
only 0.6%(//0.1), where / is the ratio of MACHO op- 
tical depth toward the L MC to that e xpected for a full 
point-lens MACHO halo (|Gouldll2005h . However, these 
~ 5(//0.1) per year events are to compared to the ~ 1 
per year currently detected toward the MCs. 

Of course, the problem of background rejection is much 
more severe because there are ~ 10 3 times more stellar 
lcnsing events. But the quality of the lightcurves is also 
higher, not only because the sources are closer but be- 
cause the bulge fields are much more aggressively mon- 
itored than the MCs, due to the possibility of finding 
extra-solar planets. This will be even more true in the 
future, as next generation microlensing planet-search ex- 
periments (already funded) begin to come on line. 

And further in the future, there are proposals to put 
a wide-field imager in space, which could obtain the 
exquisite photometry that would enable one to distin- 
guish even extremely subtle differences between point- 
lens and compact-minihalo microlensing events. Such a 
satellite would already create a syner gy between dark- 
energy and extra-solar planet studies (|Gouldll2009( ), and 
its additional power to probe for compact minihalos 
(and so the early-universe phase transitions that gen- 
erate them) would add both dark-matter searches and 
early-universe physics to the mix. 

5. CONCLUSIONS AND DISCUSSION 

We propose a new type of dark matter MACHO that 
may form copiously in the early Universe. Upper limits 
on, or the discovery of these MACHOs is a powerful new 
tool to study cosmic phase transitions and the nature of 
dark matter. The new MACHOs are ultracompact mini- 
halos that may form from cosmological perturbations in 
the Universe that are only 10-100 times larger than scale- 
invariant perturbations from inflation that seeded the 
formation of large scale structures and galaxies. These 
objects are more likely to form than PBHs, due to the 



much larger amplitude of perturbations that are required 
to form PBHs. Ultracompact minihalos may be sites of 
formation of the first Population III stars. We show that 
the magnification light curve of ultracompact minihalos 
is similar to the one of a point mass but can be easily 
recognized from it by sampling the light curve for a suf- 
ficiently long time after the peak of the magnification 
event. 

Similarly, if PBHs exist and are only a fraction of the 
dark matter, their light curve would differ from those due 
to point masses because PBHs seed the growth of an en- 
veloping minihalo that is sufficiently massive and concen- 
trated to modify the microlensing magnification curve. 
This modification is substantial whenever the mass of 
the minihalo within the lens Einstein radius is ~ 30% of 
the PBH mass. However, the PBH mass can be much 
smaller than the minihalo and MACHO masses. 

For the mass profile of the minihalo-enveloping PBHs 
as predicted by MOR07, the mass within the Einstein 
radius is only 3% of the PBH mass, thus the PBH 
magnification curve does not show measurable modifica- 
tions. However, MOR07 assumed secondary infall from 
an expanding uniform Universe seeded by a point-mass 
PBH. The later assumption in most cases is incorrect: 
PBHs form from collapse of linear fluctuations on Hori- 
zon scales, but their mass is typically smaller than the 
Horizon mass: M PB h = fH or M Hor , with f Hor < I. 
Thus, the mass that seeds the secondary infall is larger 
than the PBH mass. This leads to more massive and 
more compact minihalos enveloping PBHs, which have 
observable signatures on the microlensing light curves. 

The mass of the lens produced by compact minihalos 
scales with the mass of the perturbation as 5m 6 / 5 . Thus, 
statistical fluctuations of the mass of the perturbations 
5m produce a mass range of the gravitational lenses (i.e., 
a range of MACHO masses). This is also the case for 
minihalos surrounding PBHs, if they exist. In this second 
case we expect variations of the ratio of MACHO masses 
to PBH masses. 

Prior to this work, microlensing searches for PBHs 
were r eaching a dead-end . Although the MACHO experi- 
ment (jAlcock et al.ll200"o| ) had initially found a microlens- 
ing signature toward the MCs that was too large to be 
produced by known populations of stars, and indeed was 
difficult to explain by any objects other than PBHs, sub- 
sequent resul ts reported by EROS (iTisserand et al.ll2007f) 
and OGLE ()Wvrzvkowski et alj I2009D have tended to 
cast doubt on MACHO 's claim. But even if these are 
ultimately verified, ROM08 have argued that the density 
of PBHs required to explain the MACHO results would 
give rise to other signatures that have not been observed. 
This seemingly led to limits of order 1% on PBHs, far 
below the reported MACHO value, but also (and more 
importantly), well below the background level expected 
from Galactic and MC stellar microlenses. Hence, it had 
appeared as though PBHs were inaccessible to microlens- 
ing experiments and, in any event, cosmologically unim- 
portant (even if perhaps interesting). 

However, we have shown that PBH-seeded minihalos 
remain viable candidates for the origin of the signal re- 
ported by MACHO, and further, that even if this re- 
ported signal proves spurious, it is still possible to search 
for PBHs in microlensing data at lower levels (well be- 



low the previously conceived "floor" of the stellar back- 
ground). Both these realizations derive from the fact 
that PBHs make up only a small part of the mass of 
PBH-seeded minihalos. Hence, the cosmological density 
of PBHs can be low while the halos they seed can gener- 
ate substantial microlensing signal. And the fact that 
minihalo+PBHs are not point-like leads to a qualita- 
tively different microlensing signal. 

Moreover, we have shown that the same type of early- 
universe perturbations that give rise to PBHs, will (at 
lower amplitude) give rise to minihalos that lack PBHs. 
Such minihalos are both more likely (because they are 
produced by a larger range of initial conditions) and eas- 
ier to detect (because they deviate more strongly from 
point-lens profiles). Therefore microlensing experiments 
(past and future) are a far more powerful probe of early- 
universe physics than previously understood. 

Probing the dumpiness of dark matter at mass scales 
accessible to microlensing experiments has several impli- 
cations to understand the nature of dark matter and the 
high-energy Universe: 

1. It provides the best test to determine whether the 
dark matter is cold or warm. The masses of ther- 
mal or non-thermal relics and their free-streaming 
length could be constrained much better than in 
studies based on the power spectrum at small scales 
from the Lyman-a forest and the number of satel- 
lites in the Milky Way. 

2. It provides information on the physics of the high- 
energy Universe by constraining the amplitude of 
inhomogeneities created at phase transitions due 
to bubble nucleation, topological defects formation 
and other non-Gaussian processes. 

3. The dumpiness of dark matter in the Milky Way 
halo affects the flux of dark matter particles on 
earth, important for detections by ground based 
experiments. Unless a clump happens to inter- 
sect the earth, the average dark matter density on 
earth would be reduced with respect to the stan- 
dard value ~ 0.008 M Q pc -3 . 

4. The centers of ultracompact minihalos may be 
powerful emitters of gamma rays due to self- 
annihilation of WIMPS. A test of our model would 
consist in follow up observations of gamma rays to- 
ward the line of site of the magnification event. 

5. Ultracompact minihalos may have been the sites of 
formation of the first stars formed in our Universe. 

A positive identification of PBHs as MACHOs would 
motivate future space missions to detect signatures of 
the energy injection by PBHs in the early Universe 
(ROM08). The best measurements of the CMB spec- 
trum to date are by FIRAS on board of COBE, now 15 
years old. A follow up mission would be of great scientific 
value and easily justified if evidence is found of the exis- 
tence of non-baryonic MACHOs. Early energy injection 
also modifies the cosmic recombination history and thus 
the spectrum of anisotropics of the CMB. The Planck 
mission will improve existing upper limits on PBHs and 
may detect signatures of their existence. 
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